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Abstract 

to  the  first  part  of  the  paper,  we  discuss  existence 
results  for  a  General  stochastic  Navier-Stokes  equation  (SNS)  derived 
the  tat  Principles.  In  the  second  part,  we  deal  with  the  prop- 
" opener  Los  by  the  SNS  and  its  relation  to  statical 

moments  of  the  solution. 

Key  words:  Stochastic  Navier-Stokes,  turbulence,  Wiener  chaos,  mo- 
ments. 


1  Introduction 


Relation  of  the  Navier-Stokes  equation  to  the  phenomenon  of 

ZXSLZ  ground 

driven  by  white  noise  type  random  orce.  .  ,  .  ,  r^i  rgi  rgj 

tially  developed  and  extended  by  many  authors  (see,  e.g.  [2],  [3J,  p],  l  J. 

Pi 'to'ttee'ple^'some  form  of  stochastic  Navier-Stokes  (SNS)  equation 
was  postulated  at  the  inception  p*tA '  Xt  the 

"  the  stochastic  dimorphism 

I-Th?  wort  was PparttalIySUsupported  ^by  NSF  Grant  DMS-98-02423,  ONR  Grant 
N000 14-97- 1-0229,  and  ARO  Grant  DAAG55-98-1-0418. 
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fi(t,x)  =  n(t,v(t,x))  +  <7(t,ri(t,x))°W,  >)(0,i)=i  (ll) 

with  undetermined  (random)  local  characteristics  u(t,l)  and  a  (t,x)  .  In 
T(t,  x)  o  W  models  the  turbulent  part  of  the  velocity  field  whde 
Tmtdeh  ts  regular  component.  Following  the  classics  scheme  of  the 
"  wtonr^  ftuid  mechanics  (i.e.  coupling  (1.1)  with  Newton's  second  law) 
a  general  SNS  equation  (see  (2.3)  below),  was  derived.  It  includes  as 
sDecial  cases  the  classical  deterministic  Navier-Stokes  and  Euler  equation  as 
well  as  most  of  the  variations  of  the  SNS  equation  considered  in  the  literature. 

iTthe  present  paper  we  will  discus  new  results  on  existence  and  umque- 
ness  of  local  and  global  (pathwise)  solutions  to  these  equations i  in  the  Besse 
classes  #71,  In  addition,  we  will  derive  a  deterministic  parabohc  system  for 
the  Hernfite-Fourier  coefficients  in  Wiener  chaos  expansion  of  u  (t,  x)  which 
t  referTo  as  -  propagator” .  It  will  be  shown  that  the  statistical  moments  of 
the  velocity  field  u  (fc.x)  can  be  expressed  straightforwardly  via  the  solution 
If  the  propagator,  mile  still  an  infinite-dimensional  system,  the  Pr°P^t°r 
f  tL  qvs  eauation  is  a  much  more  simple  object  then  the  related  Kol- 
^otv  ^uaSn  On  tha  other  hand,  it  is  quite  sufficient  for  dealing  with 
basic  statistical  properties  of  solutions  to  the  SlsS  equation. 

2  Stochastic  Navier-Stokes  and  Euler  Equa¬ 
tions 

Let  (ft  T  P)  be  a  complete  probability  space  and  Y  be  a  separable 
sP-  The  scalar  product  in  V  will  *  denied  by  x  ^be  «£ 
valued  cyhndrical  Browman  motion  on  (ft,  ^  )- 

c— algebra  f I  fr(W(s),  s  <t  +  e)  will  be  denoted  >{  • 

Let  us  assume^  that  the  stochastic  fluid  flow  is  given  by  equation  (1.1) 

where  a\t,x)  and  u\t,x)  are,  respectiveiyan  V-vaiuedanda  re  -  ue 
^-adapted,  functions  on  ft  x  [0,oo)  x  Rd  .  Throughout  what  foHows the 
symbol  odW  (respectively,  -dW)  indicates  that  the  integral  is  understood 
in  the  Stratonovich  (respectively,  Ito)  sense.  F  stochastic  in  egr  ^ 
derstood  in  the  sense  of  [9],  [10],  (The  Kunita  and  Watih .m 
as  the  integrals  with  respect  to  spatially-homogeneous  Browman  motion  ar 
included  as  special  cases  of  the  above  setting). 
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In  contrast  to  the  classical  assumption  that  the  velocity  field  u  (t,  x) is  dif¬ 
ferentiable  in  t,  we  wifi  suppose  that  it  is  a  semimartingale.  More  specifically, 

we  will  assume  the  following:  w 

HI.  dn  (t,  x)  «  a  0 1 ,  x)  dt +  p  (t,  x)  •  W(t)  where  a  and  /?  are  . ^  - 

adapted  vector-functions  on  Q  x  (O.oo)  x  &  taking  values  m  R1  and  Y  , 

^Notelhat  as  in  the  classical  theory,  the  components  of  the  velocity  field, 
a  3  a  are  not  supposed  to  be  given,  and  will  be  determined  ater  from  to 
the  conservation  of  momentum  principle.  The  only  additional  assumption 
we  make  regarding  these  functions  is  that  they  are  sufficiently  smooth  m  x 

and  the  related  integrals  are  defined.  *  , 

By  the  Newton  2nd  law,  we  have  17(f)  =  F  M(*))  where  F(t,x)  is  the 

total  force  applied  to  the  fluid  particle. 

Since  the  acceleration  and  the  force  must  have  the  same  structure  we 
conclude  (for  more  detail,  see  [11])  that  there  exist  Tt  ap  e  unc  10ns 
F0  ,  F9  and  F£,  so  that 

f  tp  ({)  F  (£,  x)  dt  =  f  <p  (t)  (F0  (t,  x)  dt  +  Fa  (i.  x)  •  dW(t))  ^ 

-  f  (t)  (Ft  (t,  x )  ■  dW(t)  +  (1/2)  diFt  (i,  x )  •  F't  (t,  x))  dt. 

Similarly  to  the  classical  setting,  we  assume  that  the  forces  acting  on 
the  fluid  particle  include  pressure  and  body  forces  More  specifically  we 

assume  that  F0(i,x)  =  -Vp.(t,x)  +  €»(«,*),  G(t,x)  =  “Vp(t,x)  r  jU 
and  D (t  x)  =  -Vpt(t>x)  +  ft(i,x).  The  components  of  the  body  force  are 

considered  to  be  given  while  the  components  of  the  pressure  are  subject  o 

detI^w^tshown  in  [11]  that  (1.1),  (2.2),  and  Newton’s  2nd  law  yield  the 
following  equations  for  the  components  of  the  velocity  fie 

(  dtu  =  dj{cF& u)  -  ujdjU  -  (Va4)  •  dip  +  f  -  Vp+  ^ 

\  (g  -  Vp  -  a'diVt)  -  W,  div  u  =  0,  u(0,  x)  =  uo(x) 


a  ( t ,  x)  =  —  Vp(t,  x)  -1-  d(t,  x) 

where  p  =  p.  -  <r‘  •  ftp.  f  -  t,  -  <H  •  ftg,  and  <■«  =  K  ^  “  tie  case  °f 
ideal  Arid  L  ««  -  4  +  J.*  •  ^  if  the  fluid  is  viscous  with  the  vuscos.ty 

constant  i/. 
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3  Existence  and  Uniqueness  of  Solutions 

In  this  section  we  will  discuss  the  solvability  problem  for  the  equation  (2.3) 
in  the  Bessel  classes.  The  following  notation  will  be  used  in  the  future: 

Let  p  G  (2,oo)  and  n  €  (—00, 00). 

h;  =h;  (Rd)  is  the  space  of  generalized  functions  u  so  that  |u|„tP  = 
|(1  -  A)n/2it|P  <  oo,  where  |  •  |p  is  the  Lp-norm; 

Hn  =  H"(Rd)  is  the  space  of  all  vector  fields  u  =  ,ud)  such  that 

(E>X)l/p  <  °°- 

We  will  omit  the  sub-  and  super-indices  n  if  n  =  0. 

L„(y)  is  the  space  of  vector  functions  with  F-valued  components  gl  such 

that  MP  =  (Et  \9l\lf,p  < 

Lp  =  Up, 

55°(Rd)  is  the  space  of  all  vector  fields  <f>  =  ,4>  )  such  that  G 

Cq°  (Rd)  and  div  <f>=  0; 

The  scalar  product  in  L2  is  denoted  (•,  -)0 . 

We  will  need  the  following  assumptions: 

Bl. There  exist  constants  0  <  K  <  00  and  S  >  0  so  that  P-a.s. 

2 

(|da*7(t,  x)|  +  \aij(t,x)\  +  £  |dV(t,  x)|y)  <  K, 

k—0 


and 

K\X\2  >  [aij(t,x)  -  ^a*(t,x)  •  ai(t,x)]XiXi  >  <5|A|2 
for  all  t  >  0  and  x,  X  G  Rd. 

B2(p).  For  each  T  >  0  there  sure  constants  (7i  and  C2  such  that  P-a.s. 
for  each  s  <T  and  x  G  Rd, 

2  2 

[ dz9(s,x)\y  <  C\(s,  x)  and  ^|5*f(s,x)|  <  C2(s,x ) 

fc=0  k=0 

a^d  |Ci(t)|p  +  \C2(t)\*  <  K. 

The  following  two  results  establish  strong  (pathwise)  existence  and  unique¬ 
ness  of  local  and  global  maximal  solutions  for  the  stochastic  Navier-Stokes 
equation  (2.3) . 
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Theorem  3.1  a)  Let  p>  d.  Assume  Bl,  B2(p)  ,  and  suppose  that  |u 
oo  P-o.5.  Then  there  is  a  unique  predictable  sopping  time  C  P(C  - 
1,  such  that  for  each  stopping  time  5,  (0,  S]  C  [0,  C)  */  and  onlV  'f  ^r< 
E\-valued  continuous  solution  to  (2.3)  on  [0,5];  Also,  there  is  an  Mp- 
continuous  process  u(f)  on  [0,0  such  that  u(C~ )ii,p  =  oo  on  {C  <  oc 
for  each  S  such  that  [0, 5]  C  [0,0,  u(i  AS)  is  a  unique  solution  to  (Z 
[0,5]  •  Moreover,  for  r/=curlu,  P-o.5.  we  have 

[S\Vr!(s)\*ds<co. 

Jo 

The  function  u(t)  and  the  stopping  time  C  are  usually  referred  t 
maximal  solution  and  its  explosion  time. 

The  proof  of  the  Theorem  is  quite  involved  and  could  not  be  pre 
in  this  paper.  The  interested  reader  could  find  a  complete  proof 
forthcoming  paper  [12]. 

Now  let  us  specialize  to  the  2  D  case. 

Theorem  3.2  Suppose  that  the  assumptions  of  the  Theorem  3.1  ho 
addition,  let  us  assume  that  d  —  2,  B2( 2)  holds,  and  P-a.s.  |uo|i,p  +  |u 
oo  P-a.s.  Then  there  exists  a  unique  continuous  HpnH \-valued  solutic, 
to  (2.3)  on  [0,  oo) .  Moreover,  for  each  T  >  0 

Esup(lu(t)|ip  +  |u(t)|ij2)  <  oo 


provided  that 

E(|uo|i  p  +  luofo)  <  °°- 

Proof.The  existence  of  a  maximal  solution  follows  from  Theorei 
It  remains  to  prove  that  P(C  —  oo)  =  1.  Let  yt  =  |u(f)|i>p  +  |u 
Rm  =  inf(t :  yt  >  m)  AC-  Since  the  sequence  Rm  ’’announces”  the  predi 

stopping  time  £,  for  each  T  >  0 

P  (Rm  <T)<  P{yRm  A  T  >  m)  <  m^Ey^r, 


It  can  be  shown  (see  [12])  that 

supEyj^AT  <  oc. 
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So  limmP (Rm  <  T)  =  0,  and  therefore  P(£  =  oo)  =  1.  For  the  moments 
estimates  see  [12]. 

Remark.  One  can  consider  a  more  general  equation 

dtul  (£,  x)  =  di  (a'J  ( t ,  x)  djU1  (£,  x))  —  uk  (t,  x)  dku\t,  x) 

—dip  (t,  x)  +  bl{t,  x)diu(t,  x)  +  dip(t,  x)hlii(t,  x) 

+fl  (£,  x,  u  (£,  x)  ,  Vu (£,  x))  + 

<  (3.4) 

+[cr‘(£,  (£,  x)  +  y  (i,  x,  u  (£,  x))  -  (£,  x)]  Wt, 

divu  =  0, 

.  u(0,x)  =  u0(x),  l  =  1 

where  the  free  forces  axe  functionals  of  the  solution.  The  results  similar  to  the 
above  hold  for  equation  (3.4)  .Of  course  in  this  case  additional  assumptions 
on  regularity  of  fl  and  gl  with  respect  to  u,  Vu  as  well  as  the  appropriate 
growth  conditions  must  be  introduced  (  see  [12].  ) 


4  Wiener  Chaos  Expansion  and  Statistical  Mo 
ments 

In  this  Section  we  investigate  how  the  SNS  equation  (2.3)  propagates  chaos 
generated  by  the  driving  Brownian  motion.  Then  we  will  apply  the  Wiener 
Chaos  expansion  to  derive  formulas  for  the  statistical  moments  of  the  velocity 
field  u.  For  the  sake  of  simplicity,  everywhere  in  this  section  it  will  be  assumed 
that  the  cr,  f ,  g,  and  uo  are  nonrandom  and  the  assumptions  of  Theorem  3.2 
are  in  force. 

To  begin  with,  we  shall  introduce  the  Wiener  chaos  generated  by  W. 

Let  us  fix  a  positive  number  T  <  oo.  Let  {mk,'k  <  1}  be  an  orthonormal 
basis  in  L2  (0 ,T)  and  {£*,  k  >  1}  an  orthonormal  basis  in  Y.  Write  Mk  (i)  = 

Jo  Txii  (s)  dwk  ( t )  where  wk  ( t )  =  (W  (f)  ,  ik)Y .  Below  in  this  section  we  will 
also  use  the  notation  =  Mk  ( T )  .  Let  a  =  {af ,  k  =  0, 1, 2, ...;  i  —  1,2, ...}  , 
be  a  multiindex,  i.e.  for  every  (i,  k) ,  ak  €  JV  =  (0, 1, 2, . . . }.  We  shall 
consider  only  such  a  that  jor|  =  ^2k  i  a*  <  oc,  i.e.,  only  a  finite  number  of  ak 
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is  non-  zero,  and  we  denote  by  J  the  set  of  aJl  such  multiindice 
write  Cq  :=  TT”o  Hak  (£*)  where  Hn  is  the  nth  Hermit 

The  random  variable  =  Ga/y/ai  is  often  referred  to  as  0th 
(polynomial ). 

Let  (ej,  i  €  N}  be  an  orthonormal  basis  in  L2.  Since,  by  I 
Martin  Theorem  (see,  e.p.[13]) ,  {£Q,  a  €  J)  is  an  orthonormal 
L2  (Cl,  T?  ,  P)  ,  we  have  that  {ej  €  N,a  €  J}  is  an  orth< 

in  Li  (Cl,  Tt,  P;L2) .  This  implies  in  particular  that  for  every  v  € 
we  have  the  following  Wiener  chaos  expansion  : 

where  vQ  =  £E[vCq]  =  £  E[(v,  ej)Ca]ej.  We  will  refer 

vmnormalized  Hermite-Fourier  coefficient  of  v  (with  respect 
{ej  @€Qlz  €  N,ac  €  J})  or  simply,  Hermite-Fourier  coefficient. 

By  Theorem  3.2,  Esup3<r  |u(£)|f>2  <  00.  Thus  the  solution  c 
the  Wiener  Chaos  expansion  u  (£,  x)  =  (*■  x)  C a-  Of  cou 

problem  of  interest  is  how  to  characterize  the  Hermite-Fourie 
uQ  (£,  x) .  It  will  be  shown  below  that  these  coefficients  verify  a 
linear  parabolic  system  of  equations.  To  formulate  this  stateme 
cisely,  we  need  some  more  notation. 

Write 

(u  (t)  *  dju)Q  = 

-  Epzj  E0</J<a  ii  (Pp0)  (P+p-0)P%+a-0  (t)  9iUp+0  (t) 

Mua(t)  =  -aj(t)djUa  ( t )  -f  /{a=0}<7  (*) , 

For  v  GL2(Y),  write  £7J(v)  =  dj  f  TXi(x—y)vi(y)  dy  where  T  (x  — 
J/I/2-7T.  The  operator  Q  =  is  often  referred  to  as  gradiej 

on  L2(y).  It  is  well  known  that 

L2(y)  =  ^(L2(y))©«s(L2(y)) 

where  <S(L2(y))  =  {g  €  L2(y)  :divg  =  0}. 

For  a  €  we  define  multiindex  a  (i,  j )  6  J  by  the  formula 


R.  Mikulev ictus  and  B.  Rozovskii  / Stochastic  Navier-Stokes  Equations 


i.e.  the  multiindex  a(i,j)  might  differ  from  a  only  by  its  (i,j)  entry  which 
is  equal  to  (aj  —  l)  A  0.  Finally,  write 

DM  (uQ  (£))  =  (/af/{lQ,=1}  -  c*(t)dju a(i,k)  (£))  mi  (t) . 

Theorem  4.1  Let  u  be  the  maximal  solution  of  the  stochastic  Navier-Stokes 
equation  (2.3) .  Then  the  Fourier- Hermite  coefficients  ua  are  continuous 
¥L\-valued  functions  on  [0,  oo)  and  for  each  T  >  0, 

SUPGUaWll-p-b  |ua(t)|i,2)  <<» 


Moreover  ,  (ua  (t,  x)  ,a  E  J}  is  the  unique  solution  of  the  system 

'  (Ua  (t) ,  4>) 0  =  /{a=0}  (uo,  4>)q  +  /0‘{-(aij'3fua  (s) ,  dj<p)0- 

((u  (s)  *  diVL  (s))a ,  <f>)0  +  C V^CO)  -  9\M  (ua  ( t )) ,  <f>)0+ 

<  (4-8) 

L {a=o}  (f  (s) ,  <t>) o  +  ( dm  (ua  (s)) ,  <£)0}ds,  div  ua  =  0, 
for  all  <t>  E  Sq°  (Rrf)  and  a  E  ff. 


Sketch  of  Proof.  The  first  part  of  the  statement  follows  from  Theorem  3.2. 
It  can  be  shown  that  the  relation  E  (d,u  (£)  ux  (£)  Co)  =  (u  (t)  *  djU)Q  follows 
from  the  well  known  formula 


C.Ca  =  X) 

p<aA/? 


pKa+0-2p- 


Write  Ca  (t)  =  E  [Cal^t]  •  Note  that  Cq  (£)  verifies  the  equation 

d£ q(£)  =  mi  (£)  af  Ca<«,*) (t)dwk  (£) .  (4.9) 

Now,  the  equation  (4.8)  can  be  derived  by  differentiating  the  product 
u  (f,  x)  Qa(t)  by  Ito  formula  and  taking  expectations  of  both  sides  of  the 
resulting  equation. 
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Making  use  of  the  Wiener  chaos  expansion  (4.5)for  a  solution  of  the  SNS 
(2.3),  one  can  immediately  compute  the  first  two  moments  of  the  solution 
via  the  Hermite-Fourier  coefficients  provided  by  the  equation  (4.8)  for  the 
propagator.  Indeed,  since  EC*  =  0  for  a  ^  0  and  ECo  =  1  where  0  is  the 
zero  length  element  of  S  ,  we  have 

Eu  ( t ,  x)  =  uo  (£,  x) . 

By  (4.5)  and  Parceval’s  identity,  one  has  that  for  almost  all  x,  y  €  and 

i,s€[0,Tj, 

Eu  ( t ,  x)  u  (s,  y)  =  ^2  ua  ( t ,  x)  ua  (s,  y) . 

Moreover,  we  have 

E|u(t)£2  =  £  I** ML' 
aeJ 

Similarly,  given  the  solution  of  the  equation  (4.8) ,  the  higher  order  moments 
of  the  solution  to  SNS  equation  (2.3)can  be  obtained  by  computing  the  mo¬ 
ments  of  the  Wick  polynomials  C*-  For  example, 


Eu’  (£,  x)  v?  (£,  y)  uk  (t,  z )  = 

^  _ tfc  (£,  x)uj0(t,y)  uk(t,  z) _ _ 

o  ((o  +  P  -  7)  /2)!  ((<*-/?  +  7)  /2)!  ((/?  -  a  +  7) 

where  J>  =  €  J  :  a  +  P  - 'y  =  2p,p  €  J,0  <P  <  a  A p}  . 
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